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We explore a nonlinear field model to describe the interplay between the 
ability of excitons to be Bose condensed and their interaction with other modes 
of a crystal. We apply our consideration to the long-living paraexcitons in CU2O. 
Taking into account the exciton-phonon interaction and introducing a coherent 
phonon part of the moving condensate, we solve quasi-stationary equations for 
the exciton-phonon condensate. These equations support localized solutions, 
and we discuss the conditions for the inhomogeneous condensate to appear in 
the crystal. Allowable values of the characteristic width of ballistic condensates 
are estimated in the limit T — > 0. The stability conditions of the moving conden- 
sate are analyzed by use of Landau arguments, and Landau critical parameters 
appear in the theory. It follows that, under certain conditions, the condensate 
can move through the crystal as a stable droplet. To separate the coherent 
and non-coherent parts of the exciton-phonon packet, we suggest to turn off the 
phonon wind by the changes in design of the 3D crystal and boundary conditions 
for the moving droplet. 
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1 Introduction 



Nowadays, there is a lot of experimental evidence that paraexcitons in CU2O crystals can 
form a strongly correlated state, which can be assigned to the excitonic Bose Einstein 
condensate (BEC) 0,0,0. 

Surprisingly enough, a cloud of excitons that seems to contain the condensate can be 
prepared in a moving state, and such a packet moves ballistically through the crystal at 
T <T C . At T > T c , however, the excitonic packet exhibits the standard diffusive behavior. 
Thus, in the there dimensional CU2O crystals, the excitonic Bose Einstein condensate is 
found to be in a spatially inhomogeneous state with the well defined characteristic width 
L c h in the direction of motion [|l|] . The registered ballistic velocities of the coherent exciton 
packets turn out to be always less, but approximately equal to the longitudinal sound speed 
of the crystal, v < c s . Note that paraexcitons in the pure CU2O crystals have the extremely 
large lifetime, r ~ 13 /is, and, at the conditions we discuss in this work, they are optically 
inactive |4j]. Although transferring of the moving coherent excitonic field into a coherent 
photon field is not a hopeless task no convincing experimental results are obtained 

yet, EM- 

To understand the physics of anomalous excitonic transport, we accept that 
the macroscopic wave function \l/ ~ 0o e liPc can be associated with the coherent part of the 
excitonic packet. (Here tp c is the coherent phase of the condensate.) In other words, the 
experimental results [|I|] suggest the following decomposition of the density of excitons in 
the packet, 

n(x,t) = n co h(x,t) + An(x,t), (1) 
where n coh (x, t) « n coh (x — vt) is the ballistic (superfluid) part of the packet, 

n coh (x-vt) ~\y \ 2 (x-vt), (2) 

and An(x, t) is the diffusive (non-condensed) part of it, 

An(x,t) ~ (5$5ip(x,t)}. (3) 

Therefore, the problem is how to describe a spatially inhomogeneous state of the excitonic 
BEC in terms of ^(x^) and 8^>(x.,t), where 5ip is the "fluctuating" part of the exciton 
Bose field. Indeed, if a coherent excitonic packet moves in a crystal (or another semicon- 
ductor structure), it interacts with phonons, non-condensed excitons, impurities and other 
imperfections of the lattice, etc., (see Ref. [[Top). 

All this makes the problem of superfluidity of the Bose-condensed excitons a rather 
complicated and challenging one. 

In this Letter, we start from the simplest possible approximation: we describe the 



moving condensate only and show that a sort of Gross Pitaevskii equation pT]|,||12| do has 
sense at T <C T c . In this case, the finite characteristic length of the condensate (but not 
of the total exciton-phonon packet) appears naturally in the framework of the effective ID 
nonlinear Schrodinger equation, by which we model the real 3D conditions. 
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2 Exciton-Phonon Condensate 



To obtain the necessary density of excitons in the excitonic cloud and, thus, meet the 
BEC conditions, the CU2O crystals were irradiated by laser pulses with Tiujl 3> E gSuP at 
T ~ 2 ~ 5K. Note that the cross-section area of an excitation spot on a surface of the 
crystal, S, can be made large enough, so that S ~ S sm {. Then, the classical "phonon wind", 
or the flow of nonequilibrium phonons from the surface into the bulk JT3, can transfer the 



nonzero momentum to the excitonic cloud, P cxc ~ iV x (h~k ) 7^ and P cxc _L S Bm f. 

As a result, the packet of moving excitons and nonequilibrium phonons of the phonon 
wind (iVph ~ N x ) is actually the system that undergoes the transition toward the Bose 
Einstein condensation. Then, one can estimate the energy of the packet without any con- 
densate as follows 

E ~ N x {(h 2 /2m x )(k 2 Q ) + 2u n x } + N ph (WJ, 

where h(ko x ) = m x (v) ~ m x c s and Uq > is the exciton-exciton interaction strength, and 
n x is the average density of excitons in the packet. 

At T < T c (n x ), or, equivalently, n x > n c (T) Q, the condensate can be formed inside 
the excitonic droplet and the following representation of the exciton Bose-field holds, ip = 
\l/o + <H^- (As the kinetics of condensate formation is not a subject of this paper, we 
assume T cloud ~ T and T — > 0.) Moreover, for the displacement field of the crystal, we can 
introduce a nontrivial coherent part too, i.e. u = u + Su and u 7^ 0. Then, a moving 
coherent packet can contain the exciton-phonon condensate, which is the self-consistent 
exciton-phonon field, \l/o( x >£) • u o( x >^)- 

The macroscopic wave function of excitons, \l/o( x , t) ~ ^0(^5 1), is normalized as follows 

[\% \ 2 (x,t)dx = sj<j>l{x)dx = N , (4) 

where N Q is the (macroscopic) number of condensed excitons, and, generally, N Q 7^ N x . To 
model the ballistic motion of n co ^(x,t), we use the following ansatz, 

*o(x,t) = e-^ + ^l 2 -^ 1 S^ +k ^ <p o ( x - vt), (5) 

u j{x, t) = u (x - vt) 8ij, (6) 

where E g = E gap — E x , E x is the exciton Rydberg, if = const, hk = m x v, and /i is the 
effective chemical potential of the condensate. Note that there is no difference between the 
ballistic velocity of the exciton-phonon packet and the superfluid velocity of the condensate 
if we start from ansatz (^). 

For the envelope function <f) (x) and the phonon direct current u Q (x), we obtain the 
following stationary equations [|] . || 

- \fM\Mx) = (-(ft 2 /2m x R 2 + u (j>l(x) + ^^(x))0 o (x), (7) 
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d x u a (x) ~ consti <j) Q {x) + const2 4>o{x). (8) 

Here, the "bare" exciton vertices, both the two-particle v > and the three-particle 
v\ > 0, can be strongly renormalized because of exciton-phonon interaction. We choose it 
in the simplest form of Deformation Potential, 

-> {E g + aodju^m, a > 0. (9) 

Yet the first (cubic) anharmonicity K3 7^ is taken into account to model the lattice. 
Note that we consider the case v ~ c s , and, at first glance, the (adiabatic) assumption of 
coherent propagation of the crystal deformation field and the excitonic condensate might 
not be valid. However, the nonlinear lattices support some localized excitations that can 
consist of two parts, namely, the direct current and alternating one, 

u ~ u Q (x — vt) + 5u(x — vt, t). 



Such excitations can move with the (group) velocity v ~ c s ||14| . Therefore, Eq. (|8|) can be 
used to describe the exciton-phonon condensate with v ~ c s . 
If the following conditions 

z>o = ^0(^0, &o, v/cs) < and u t = v x {v x , a , v/c s , « 3 ) > (10) 

can be valid, the localized (solitonic) solution of Eq. (0) exists. It can be written in the 
following form: 

<j> (x) = $ f(J3($ )x, const($ ) ), /3($ ) = V /(2m x /^ 2 )| yU |($ ), (11) 

where $ Q is the amplitude of the soliton and 1/ f3(§ ) = L is its width. Then, the charac- 
teristic length of the condensate can be estimated as 



L ch ~ (2 ~ 4) L oc 



-1/2 



and the effective chemical potential of the condensate, \fj,\(N , v), defines the value of L c h 
in this model. 

We can use the simplest approximation, 

H($ )«|£o|$o/2 and \n\(N , v) oc iV D 2 , (12) 

which turns out to be valid at n a > 5 — 10. Here n Q = N*/N Q , and N* = 2S/a% is a 
macroscopically large parameter (a x is the exciton Bohr radius), 

N* ~ 10 13 ~ 10 14 . 

Indeed, one can estimate more accurately, 

ImI(*o) = l^ol *o/ 2 - ^1 *o/3 and \^\(N ,v)(x(2n 2 xE x al + Aiy^v/c^y 1 . (13) 
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Here, x = /i x /m x and the last estimate is valid for iV < N*. 

We choose approximation (|12|), and the characteristic length of the condensate can be 
estimated as follows ( u$ = Eo(v/c s ) a x and z/ = Sq a x ): 



L ch (N , v) 



Ti 2 . . E. 



$- 2 (iV , v ) ~ 4 — * n Q a x oc iV- 1 . (14) 



At u = f Q , or, equivalently, at 7(f) = j , where 

j(v) = c 2 /(c 2 — v 2 ) and v < c s , 

we have eo(v Q /c s ) = and 7 Q ~ 3 ~ 5 |§. Therefore, solitonic solution (|ll|) and estimate 
( |T4| ) are valid at t> < t> < c s , or 7(1;) > 7 Q . For example, at 7(f) ~ 2 7 Q ~ 6 ~ 10, 
(7(1* = 0.95 c s ) ~ 10), we obtain i> < and estimate Eo(v/c s ) ~ — £o- 

Thus, at 7(f) > 7 and n Q > 10 (for estimates we take \so(v/cs)\ ~ (10~ 2 ~ 10 _1 ) i? x 
and n D ~ 10 1 ~ 10 2 , i.e. iV ~ 10 11 ~ 10 12 ), we obtain a large factor multiplied by a x as an 
estimate of L ch ||, e.g., 

L c h = J 7 a x ^> a CI . 

Here a cr ~ 4 A is the lattice constant and a 3 . ~ (3 ~ 4) a 3 r for the paraexcitons in Cu 2 0. 

Within the quasi-stationary approximation (|)-(0) at T = 0, this is a physically rea- 
sonable result, T ~ 10 2 ~ 10 4 , and the duration of the condensate can be estimated as 
t ch ~ 2 ■ (HT 11 - 10~ 9 ) s. Although these results are in a qualitative agreement with the av- 
erage 3D densities of the Bose-condensed excitons in moving packets at T 7^ 0, n Q ~ 10 17 ~ 
10 18 cm -3 IJ, the characteristic duration of such localized packets is At ~ (4 ~ 8) • 10 _7 s 
experimentally. 



3 Stability against Outside Excitations 

Analyzing experimental results jj], one can notice that a long-lasting "tail" of excitons is 
followed by the coherent (localized) excitonic packet. This tail might be explained by insta- 
bility of the exciton-phonon condensate moving in the lattice. Indeed, such an instability 
can lead to continuous emission of the excitons out from the condensate. As dtN Q < 0, 
we have dtL c h(t) > and, as a result, nonstationary transport of the ballistic condensate. 
Alternatively, the condensate is stable, and diffusive propagation of the non-condensed ex- 
citons, AiV ~ iV x — N a , is responsible for the tail in n(x,t). We argue that the second 
scenario seems to be true. 

For outside collective excitations (see Fig. 1), the asymptotics of the low-lying energy 
spectrum can be found easily. Indeed, if we assume that <p 2 (x) ~ and d x u Q (x) ~ in the 
outside packet area, the excitonic and phonon branches are (formally) uncoupled. Then, 
for the excitonic branch, we obtain in the co-moving frame 

^e X (k) « + (h 2 /2m)k 2 , u k (x) « w k e ikx , v k (x) « 0, |x| > L , (15) 
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where Uk and Vk are Bogoliubov-deGennes amplitudes, and co> p h(k) = c s |k| in the laboratory 
frame of reference. Note that the condition hu cx (k) + hk x v > can be violated at the 
velocities close to v Q . This is the hint that the instability regime can occur ||. 

We assume that the exciton and phonon can be emitted from the condensate coherently. 
Then the emission of SN excitons out of the coherent packet can be described as an appear- 
ance of SN outside collective excitations in the system condensate plus medium, see Fig. 
1. Therefore, one can use Landau arguments to analyze the stability of the condensate. It 
turns out that the moving condensate is stable against the direct emission of outside exci- 
tations. However, this is not the case for the inside excitations. In the laboratory frame, 
these excitations can be represented in the form 

«(^o+kx) -i(ui k +k x v)t , / _ .\ -i(tp +kx) i(uj k +k x v)t 



<ty(x, t) ~ u k (x - vt) e ^ 0+k *> e - l[UJk+k * v)t + v k (x - vt) e"^ 0+kx) e 



and 

5u(x, t) ~ C k {x — vt) exp(i({p + kx)) exp(— i{oj k + k x v) t) + c.c, (16) 

and the following asymptotic behavior is valid for the Bogoliubov-deGennes amplitudes 
(<p o (x/L ) ~ exp(-|x|/L ) ): 

Uft(ic) ~ v fc (a?) ~ (f) o (x/L ), C k (x) ~ (f>l(x/L ) at \x\ > L . 

We can start from the standard criterion, 

H^k - \kx\v) > at l/c^ ~ zL \ (17) 

where the superscript (— ) means k x < and z ~ 3 ~ 10 corresponds to the low-lying inside 
excitations, and find under which conditions inequality (pj]) is valid. 

It turns out the analog of Landau criterion for Eq. QT7| ) is rather simple 

\p\{N , v) > /x cr ~ const (2m x c s 2 ), (18) 

where m x ~ 2.7 m e , c s ~ 4.5 ■ 10 5 cm/s. For z ~ 3 ~ 10, we estimate const(;?) in Eq. flTjp 
as of the order of 10 _1 , and we obtain the following result 

fi a ~ 10~ 5 - 10~ 4 eV ~ (10- 4 - 10" 3 ) E x . (19) 



Qualitatively, inequality flT8"| ) is valid for the relatively high velocities and numbers of the 
condensed particles. In theory, one can fix the number n Q = N*/N Q and obtain an analog 
of Landau critical velocity, but v cr < v < c s is the stability criterion in this case. 

However, within approximation (12"D,fll4D and with |e (^/c s )| — 10 _1 E x and n Q ~ 10 1 , 



we have |/i| ~ /i cr and L c h — 4 • 10 2 a x ~ (2 — 3) • 10 3 A. We speculate that such a ballistic 
condensate can be considered as a stable one in the limit of T — ► 0, see Fig. 1. For 
comparison, the condensates with |e | ~ lO -2 ^ and n Q > 50 — 100 seems to be unstable 



6 




1-outside-exc *1 < *2 v > v t 



Figure 1: The ballistic condensate, 4> (x — vt) ■ u (x — vt)8\j, seems to be stable in relation 
to the direct emission of outside exciton-phonon excitations. (We consider the backward 
emission in the long-wavelength limit, \k x \ ~ z Lq 1 and z <C 1.) The outside excitations 
presented on this figure are labeled by the wave vectors, k x , k' x < in the co- moving frame. 
To a first approximation, the outside excitations can be described in terms of free excitons 
and free (acoustic) phonons emitted from the condensate coherently. 

against the inside excitations. If this is the case, the continuous generation of the low-lying 
inside excitation takes place, 

N ^N - 6N(t), d t 6N in . cx (t) > 0, 5iV in _ cx (t) > y/F Q . 

This process can be accompanied by the emission of outside excitations as well. 

To conclude, more accurate investigation of the stability problem is necessary. Indeed, 
the characteristic width of the quasi-stationary solution near the threshold of stability is 
of (1 — 3) • 10 3 A, whereas the typical length of the crystal used for experiments is of 
(2 — 4) ■ 10 _1 cm. We speculate that the ballistic (superfluid) propagation of the condensate 
is more than changing x — > x — v t and adding exp(ik x) to the solution of Eqs. (|7|) and 
(H). For example, the localized solitonic solution fl5|),(fj) can be used as a reasonable initial 
condition for modeling of how the coherent part of exciton-phonon packet actually moves 
in the presence of thermal phonons and/or point scattering centers, etc.. 
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4 Discussion 

Recall that the self-consistent exciton-phonon condensate seems to be only a part of the real 
moving packet. The noncoherent part of it, the noncondensed excitons An(x,t) and the 
unidirectional phonon wind Au(x,t), effects the propagation of the condensate. Here we 
address the question on whether it is possible to diminish (ideally, to turn off) the phonon 
wind after the moving exciton-phonon condensate has been formed. As a result, the diffusive 
noncondensed excitons can be delayed, and the coherent signal and the noncoherent one 
are separated in time. 

We consider a system consisting of a crystal (semiconductor) and two ideal conductors 
(metals). The geometry of such a system is shown on Fig. g. Let I = |db| = |ac| be the 
distance between two conducting planes surrounding part A'B of the crystal. Let A be a 
wavelength of the light emitted during the exciton recombination. (Usually, A lies in the 
visible wavelength area, ~ 500 nm.) If the following relation 



can be realized by appropriate changes in the growing technique, an interesting physics of 
penetration into a channel appears. Here D ~ 1 — 2 mm is the diameter of an excitation 



Note that due to boundary conditions on the conducting walls, it is impossible to emit 
light with the wave vector less than n/l and the quantum energy less than heir /I. If the last 
value is larger than the energy of the electron-hole recombination, the later is suppressed in 
the region abed. Then the lifetime of excitons (both free and bound ones) can be increased 
by several orders of magnitude. 

We are interested in evolution of the exciton-phonon system when the cloud reaches x = 
B (see Fig. 2), and the excitons begin to move in channel abed. Roughly speaking, we have 
such a connection problem, in which the boundary conditions and the characteristic length 
scales (L y in our case) can be quite different for the 3D crystal and the channel, respectively. 
As boundaries ab and cd of the channel are metallized, the boundary conditions of the 
displacement field u(x,y,t) ~ (u x ,0,0) can be taken zero along them. 

We consider transition of the phonon wind from part AB of the crystal to the channel 
by posing of the time-dependent boundary condition at the connection B | bd (see Fig. 2): 



where f(u>) is localized inside \u\ < u* . For the displacement field inside the channel, we 
can write the following representation, 



I < A < D 





(20) 




n=odd 



(21) 



S 



Obviously, all the modes in Eq. fl2"T|) with frequencies less than u ~ c s 7r// are dumped 
inside the channel. Let the following inequality be valid 

2tt 2irc s nc s 
r A / 

Here, r and A can be estimated from the duration and the characteristic width of the (total) 
exciton-phonon packet driven by the phonon wind, r (1 ~ 5) ■ 10 _7 s and A ~ 10~ 2 cm. 
Then, only the "tails" of f{w) are actually transfered into the amplitudes A n (uj) of the 
moving packet u n (x,t). Hence, the phonon drag force (oc d x u(x,y,t)) can be essentially 
suppressed inside the channel, and the non-condensed excitons can diffuse almost freely 
there. 

We speculate that the exciton-phonon condensate being formed in AB part of the 3D 
crystal can, first, penetrate into and, second, pass through the channel (see Fig. 3). 
This problem seems to be not a quasi- stationary one. Indeed, within approximation 

— i/>a( x ,y,t) an d u = (u ox (x,y,t), 0, 0), we have the following system to describe 
the dynamics at x > B: 

h 2 

i?id t ip (x,y,t) = ( - ^-A Xjy + is \il) o \ 2 (x,y,t) + z/i|^ G | 4 (x, y, t) }ip (x,y,t) + (22) 

+ (T d x u ox (x, y, t) if) (x, y, t), 



(d 2 - c 2 A Xiy )u ox (x, y, t) -c 2 2k 3 d 2 u ox djU ox (x, y,t) = p V d x (\ip \ 2 (x, y, £)), (23) 

if> o (x,-l/2,t) = ii> o (x,l/2,t) = 0, u ox (x,-l/2,t)=u ox (x,l/2,t) = 0. (24) 

Boundary conditions ([24]) model the metallized surfaces of the channel. 

Here, we model a kind of smooth penetration, for example, no vortexes are generated 
near the 3D-channel intersection. Recall that the width of the condensate in our quasi- ID 
model is L ~ (10 2 - 10 3 ) a x < A, and r ~ L /c s ~ lO" 10 - 10" 9 s< r. As L < I < A 
and uj* < ujq < t^ 1 , we assume that the localized nonlinear condensate can penetrate into 
the beginning of the channel without the loss of its coherence. Then, instead of solving 
Eq. (]22|)-(p4|) with "pumping", (i.e., time dependent) boundary conditions at x — B, one 
can transform the quasi-stationary solution (|5|) and (|5]) into the initial condition inside the 
channel and choose t — 0. For example, one can use the following trial functions: 

ip o (x,y,0) ■ u o (x,y,0)5 lj = 
= exp(im x vx) $ G cosh _1 (Lo ^ ) 4>{y) ■ [Q Q - Q D tanh(L ~ 1 x) J Q(y), (25) 
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where 4>(y) and Q(y) have to satisfy the boundary conditions (^4|) and be self-consistent 
as well. We assume that near the boundaries they have the same coherence length L 0y . It 
can be estimated as IIIB 



L l y ~ y/v \Mx)\ 2 ™x/h \J(2m x /h 2 )u Q ^j2 cos\T l {L^x ). (26) 

Near the soliton maximum (x — 0), we obtain L 0y ~ L (see Eq. ( p!4|) and z/ ~ |z> | is valid 
in the stability area). 

We estimate that as many as ~ N Q (l/D) of the iV original Bose condensed excitons 
can penetrate from the crystal into the channel if only 1/2 > L 0y (x) for \x\ < (2 ~ 4)L . 
Then, the BEC conditions can be saved for the droplet inside the channel at T ^ 0, and the 
exciton-phonon condensate can move through the channel ballistically (Fig. 3). However, 
the question of stability of such a motion remains to be open ||16|| . Experimentally, the task 
is to register the signal from the output of the channel, x = ac. For example, a small 3D 
part could be attached to the channel in order to convert the excitonic mass current into a 
measurable electric signal, see Figs. 2 and 3. 



5 Conclusion 

In conclusion, we note that the phonons play a crucial role in almost all the current models 
aimed to explain or predict coherent behavior of excitons in semiconductors, see, f.ex., 
10[ , [[HJ , ||18[ . Although our model of the exciton-phonon condensate fails to predict the 



width of the exciton-phonon packet correctly, the theory yields a qualitative description of 
the experiments and a reasonable value for the critical velocity. However, one has to take 
into account the thermal excitons (e.g., the weak tail that is always observed behind the 
soliton JIJ) and the thermal phonons of the crystal to make the model with condensate more 
realistic. 
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Figure 2: In order to separate the coherent and noncoherent components of the total exciton- 
phonon packet at T <C T c , it is crucial to turn off a phonon wind. (The phonon wind, or 
a sequence of nonequilibrium acoustic phonons, "blows" unidirectionally from surface A of 
the crystal.) This can happen when the packet n(x, t) enters the channel abed of the width 
/ in the y direction. The characteristic lengths of the total packet and the condensate, A 
and L , respectively, satisfy the inequality L < I < A. Then, the condensate could move 
ballistically through the channel with the velocity v, whereas the non-condensed excitons 
will diffuse (almost) freely inside the channel because the phonon wind cannot penetrate 
into it. As a result, the (total) excitonic current n(x,t') ~ n (x — vt') + An(x,t') from the 
output of the channel (x = ac) can be converted into the electric current near the surface 
B'; |A'B'| < |AB|. 
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Figure 3: We assume that the transition from the 3D crystal AB to the channel abed can 
be done smoothly enough in order not to create turbulence in the condensate penetrating 
the channel. Then, one can expect the solitonic shape of the coherent part of the droplet 
at the end of the channel. (The noncoherent part at t — t' is not shown on this figure.) 
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